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ABSTRACT 

It is often desirable to characterize a turbomachin- 
ery flow field with a few lumped parameters such as 
total pressure ratio or stage efficiency. Various av- 
eraging schemes may be used to compute these pa- 
rameters. This paper describes and compares the mo- 
mentum, energy, and area averaging schemes. The 
schemes were compared for two computed solutions of 
the midspan section of a transonic fan stage: a steady 
averaging-plane solution in which average rotor out- 
flow conditions were used as stator inflow conditions, 
and an unsteady rotor-stator interaction solution. The 
solutions were computed on identical grids using sim- 
ilar Navier-Stokes codes and an algebraic turbulence 
model. The unsteady solution is described, some un- 
steady flow phenomena are discussed, and the steady 
pressure distributions are compared. Despite large un- 
steady pressure fluctuations on the stator surface, the 
steady pressure distribution matched the average un- 
steady distribution almost exactly. Stator wake pro* 
files, stator loss coefficient, and stage efficiency were 
computed for the two solutions with the three aver- 
aging schemes and are compared. In general, the en- 
ergy averaging scheme gave good agreement between 
the averaging-plane solution and the time- averaged un- 
steady solution, even though certain phenomena due to 
unsteady wake migration were neglected. 

INTRODUCTION 

Numerical analyses of three-dimensional steady flows 
in isolated turbomachinery blade rows can now ac- 
curately predict blade row performance in reasonable 
computer times [1, 2, 3, 4, 5]. It is thus reasonable to 
consider approaches for calculating the performance of 
multiblade-row machines. 

One approach is to compute unsteady rotor-stator 
interaction directly, integrating the time-averaged per- 
formance from the unsteady solution. This was done 
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in three dimensions by Rai [6] with excellent results, 
but solution times were long. In more recent 3-D un- 
steady calculations, solution times have been greatly 
reduced [7, 8]. However full 3-D unsteady calculations 
still require substantial computer resources. 

Several researchers have computed unsteady rotor- 
stator interaction in two dimensions. A 2 1/2-stage 
compressor was modeled in [9], and [10, 11, 12] have all 
modeled single-stage machines. Here solution times are 
more tractable, but results still require time averaging 
for interpretation. 

Another approach is to model a multistage machine 
as steady in some average sense. Adamczyk’s average 
passage model provides a formal mathematical frame- 
work for this type of analysis and has been used very 
successfully [2]. Other researchers [3, 4, 5] have used 
a more ad hoc approach in which the solution is cir- 
cumferentially averaged at a plane between blade rows. 
This averaging-plane approach has given useful results 
in many instances, but raises the question of how to 
average the solution for the best results. 

In the present work (originally described in [13],) 
the midspan section of a transonic fan stage desig- 
nated NASA stage 67 [14, 15, 16] was analyzed using 
both a steady averaging-plane approach and an un- 
steady rotor-stator interaction approach. Three aver- 
aging schemes, referred to as the momentum, energy, 
and area averaging schemes were used to average the 
steady solution between the rotor and stator. The un- 
steady solution did not require averaging between the 
blade rows but did require averaging to show time- 
averaged quantities as a function of space (e.g., surface 
pressure distribution) or to show spatially-averaged 
quantities as a function of time (e.g., stage efficiency). 
The unsteady solution was considered a standard by 
which to evaluate the three averaging schemes used in 
the steady solutions. Comparisons were made between 
unsteady and steady pressure distributions, wake pro- 
files, stage adiabatic efficiency, and stator loss coeffi- 
cient. In general, the energy averaging scheme gave the 
best agreement between the unsteady and steady solu- 
tions. Qualitative comparisons of rotor wake migration 
were made with experimental measurements made by 
Hathaway et al. [16]. 
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FAN STAGE 

Calculations were made for the 50-percent mass-flow 
streamline of a transonic axial flow fan stage desig- 
nated NASA stage 67. The fan was designed for a 
rotational speed of 16 043 rpm and a mass flow of 34 
kg/sec. The rotor and stator are separated axially by 
about 85 percent of the rotor chord at midspan in or- 
der to minimize noise. Because of the large interblade 
spacing, the unsteady effects on the stator are mostly 
due to interaction with the rotor wakes. Although this 
allows wake-generated unsteadiness to be studied in- 
dependently of acoustic effects, it also means that the 
conclusions drawn here regarding averaging schemes do 
not necessarily extend to closely coupled blade rows. 

The midspan radius is about 19.3 cm and is nearly 
constant. The stage has a significant amount of hub 
and tip convergence, such that the annular area de- 
creases by 22 percent through the rotor and by 27 per- 
cent through the stage. The radius change and annulus 
convergence were modeled using the quasi-3-D analy- 
sis described in the following section. The rotor has 
22 multiple-circular-arc blades designed for purely ax- 
ial inflow. The stator has 34 double-circular-arc blades 
designed to produce purely axial outflow. (A controlled 
diffusion stator was also tested in [14, 16] but was not 
considered here.) 

The fan has been extensively tested using both 
conventional aerodynamic probes and laser anemom- 
etry. The ensemble average and variance of the laser 
data have been analyzed to compute the rotor-wake- 
generated and unresolved unsteadiness [14, 15, 16]. In 
the present work, qualitative comparisons were made 
to the rotor-wake-generated unsteadiness data. 

NUMERICAL METHOD 

Steady solutions were made using the RVCQ3D code 
developed by Chima [17] and unsteady solutions were 
made using the URSA code developed by Jorgenson 
and Chima [10, 11]. Both codes solve the quasi-three- 
dimensional thin-layer Navier-Stokes equations for a 
blade-to-blade stream surface and include the effects 
of rotation, radius change, and stream sheet thick- 
ness. The Baldwin-Lomax model was used for tur- 
bulent flows. The flow equations were discretized us- 
ing second-order-accurate finite differences, with third- 
order artificial viscosity added everywhere and first- 
order artificial viscosity added near shocks. The equa- 
tions were solved using an explicit four-stage Runge- 
Kutta scheme. Whenever possible, the same numerical 
parameters (i.e., artificial viscosity coefficients, turbu- 
lence model parameters, etc.) were used for both the 
steady and unsteady calculations. 

Circumferentially periodic C-type grids were gener- 
ated using the GRAPE code [18]. Grid sizes were 
303 x 45 for the rotor and 307 x 45 for the stator. 
Identical grids were used for both calculations. 

Radius change and endwall convergence were mod- 


eled in the quasi- three-dimensional analysis by assum- 
ing that the blade-to-blade flow follows an axisymmet- 
ric stream surface of known radius r and normal thick- 
ness h . The equations were formulated in an (m,0) 
coordinate system, where m is the arc length along the 
surface, dm 2 — dr 2 + cfz 2 , and r = r(m) and h = h(m ) 
must be specified. Details of the quasi- 3-D formulation 
may be found in [10, 11, 17]. 

The stream surface and radius were extracted from 
steady 3-D analyses of the rotor (see ref. 1) and the 
stator (unpublished.) The 3-D solutions were first cir- 
cumferentially averaged to produce an axisymmetric 
solution. The mass flow was then integrated as a func- 
tion of span, and the location of the 50-percent stream- 
line was tracked. Finally the blade profile, radius, and 
relative distance between neighboring streamlines were 
interpolated as functions of distance along the stream- 
line. The 50-percent streamline is near midspan where 
the radius only varies about 3 percent through the 
machine. The stream surface thickness h(m) varies 
much like the annulus area, decreasing about 27 per- 
cent through the stage. 

The steady quasi-3-D code uses a spatially- variable 
time step and implicit residual smoothing with a CFL 
number of 5.5 to accelerate convergence to a steady 
state. Stage calculations were made by solving an iso- 
lated rotor and then averaging the exit flow. The av- 
erage total pressure, total temperature, and tangential 
velocity were specified as inlet boundary conditions for 
an isolated stator solution. The stator exit static pres- 
sure was varied iteratively until the stator and rotor 
mass flows matched. Solution times were about 3 min- 
utes per blade row on a Cray Y-MP. 

The unsteady quasi-3-D code uses a constant time 
step subject to a CFL limit of about 2.8. Spatially 
varying implicit residual smoothing was used to in- 
crease the maximum CFL number to about 12 while 
retaining second-order accuracy in time [11]. The ro- 
tor and stator grids overlap by one point at an interface 
between the blade rows. The interface was updated by 
a nonconservative interpolation of the flow variables. 
Two rotor and three stator blades were solved to ap- 
proximate the 22:34 blade count of the actual stage. It 
took 1126 iterations at 0-52 seconds per iteration on the 
Cray Y-MP for two rotor blades to pass three stator 
blades. This time interval will be referred to as a blade 
passing interval . Stator lift coefficients were roughly 
periodic after about eight blade passing intervals and 
fully periodic after 16. Results are shown after 21 blade 
passing intervals (almost two rotor revolutions) which 
took 3.4 hours on the Cray. 

AVERAGING SCHEMES 

It is often desirable to characterize a complicated 
flow with a few lumped parameters. For example, a 
computational analysis of a fan stage might specify 
constant p 0 and t q at the inflow boundary and constant 
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p at the exit, and a turbomachinery designer would 
want to know the overall efficiency of the stage. In 
a real machine, each of these quantities varies in both 
space and time so that it becomes necessary to develop 
representative averages. 

In two dimensions any two independent kinematic 
properties (velocity components) and two independent 
thermodynamic properties completely characterize the 
fluid state. As such, four properties may be averaged 
to satisfy four desirable characteristics of the system 
such as global conservation of mass, momentum, and 
energy. Since flow properties are related nonlinearly, 
the average properties may not necessarily satisfy other 
characteristics that were satisfied by the original sys- 
tem; that is, information is lost through the averaging 
process. It thus becomes necessary to decide what in- 
formation must be retained for a particular application 
and to devise averaging schemes accordingly. 

Any arbitrary quantity can be averaged in a math- 
ematical sense; however, averages of some quantities 
have more physical significance than others. Flow 
properties are often categorized as being either inten- 
sive or extensive. Intensive properties are single valued 
at each point in the flow. Thermodynamic variables 
such as pressure and density are examples of intensive 
properties. Extensive properties depend on the mass 
or volume (i.e., the extent) of a system and are often 
defined as an integral of an intensive property over its 
respective mass or volume. Mass =: f pdV and force 
= f pdA are examples of extensive properties. Exten- 
sive properties have the characteristic that they may be 
summed over the parts of a system to give a total value 
of the property. This characteristic makes it desirable 
to use extensive properties for flow-field averages. 

Three averaging schemes are described below. They 
are referred to as the momentum, energy, and area av- 
eraging schemes, or MAS, EAS, and AAS, respectively. 


MOMENTUM AVERAGING SCHEME (MAS) 

The MAS is the most rigorous of the three schemes 
described here and serves as the basis for the other two. 
For 2-D adiabatic flow it is equivalent to the mixed-out 
average often referred to in turbomachinery literature. 

The flux terms in the Euler equations are intensive 
properties that may be integrated over their respective 
areas to produce extensive properties. The unsteady 
twodimensional Euler equations are given by 


where 
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and 

e + p = p(C v T + V 2 /2) + pRT = ph 0 

Assuming steady flow and integrating in y over the 
pitch / of a turbomachine gives 

^-J^Fdy + G(l)-G(0) = 0 

If the flow is periodic over /, then G(l) — G( 0) = 0, and 
I Fdy = constant (2) 

Jo 

Applying (2) to a uniform flow (denoted by an overbar) 
gives 

| J F(p,u,...)dy = F(p, «,-.) = F(p,u...) (3) 

This expression may be used to define an average flow 
by solving for p, u, . . . in terms of the local integrals 
J c , I Xy ... as follows: 
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The average quantities are given by 
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V = Iy/Ic 

(10) 

h~o = h/h 

(11) 

Using the definition of ho, 
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and eliminating p, ti, and v using (8, 9, and 10) gives 
a quadratic equation for p. The solution is 
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The positive sign is for axial subsonic flow and the 
negative sign is for axial supersonic flow. Once p is 
known, the remaining quantities may be found directly. 


Equations (4-13) define a complete set of average 
properties. Since the average flux and average proper- 
ties are constant in x (eq. 2,3), these properties rep- 
resent the uniform flow far downstream of a cascade 
after the blade wakes have mixed out, and thus the av- 
eraging scheme is the so-called mixed- out average. Far 
downstream, viscous dissipation has caused all velocity 
gradients in the wakes to decay to uniform conditions 
with a subsequent increase in entropy and decrease in 
total pressure. 

The Euler equations (1) apply to 2-D adiabatic flow, 
which may not be the case in a general turbomachinery 
problem. Nevertheless, equations (4-13) may still be 
used as an averaging scheme. The MAS is relatively 
insensitive to the axial location where the averaging 
is performed. It accurately represents available thrust 
but tends to give low efficiencies and high losses, de- 
pending on the degree of flow nonuniformity. 


ENERGY AVERAGING SCHEME (EAS) 

Total pressure and temperature are commonly mea- 
sured experimentally with combination probes located 
fairly close behind turbomachinery blade rows, well be- 
fore the wakes have mixed out. Behind a rotor, it is 
generally assumed that the measurements are somehow 
temporally averaged by the low frequency response of 
the instrumentation. This temporal average is equiv- 
alent to a circumferential average if the rotor speed is 
constant. Behind a stator, measurements are taken at 
many locations and are usually mass averaged. Total 
pressure is an intensive property and a mass average 
has questionable physical significance. However, mass 
averaging computed total pressures may give the best 
comparison with experimental data averaged the same 
way. 

In the EAS the axial momentum integral (5) is re- 
placed with a mass average of the ideal total enthalpy 
h 0 i = Po 7_1)/7 (nondimensionalized by some arbitrary 
reference state): 


Li = 



(14) 


Primitive variables are found from the integrals J c , J y , 
I t > and I c i in a manner similar to that used in the MAS, 
except that a nonlinear equation must be solved: 
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where V ? = u 2 + v 2 . 
variable: 


It is convenient to define a new 



Then p can be written in terms of $ as follows: 
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Equation (15) can be solved for $ using Newton it- 
eration. Note that the root $ must be between zero 
(Af = 0) and one (Af = oo.) 

The EAS gives a good measure of the local total 
pressure and temperature, and hence the local adia- 
batic efficiency and loss for a flow. It usually gives the 
best comparison with experimental data. Far down- 
stream in a 2-D flow the EAS gives about the same 
result as the MAS. 

In recent unpublished work, Tweedt has used an en- 
tropy averaging scheme in which the local specific en- 
tropy of the flow is mass averaged instead of the ideal 
total enthalpy. The entropy averaging scheme is consis- 
tent with the second law of thermodynamics, and thus 
gives the total pressure which best represents the loss 
at a given axial location. In practice, the energy and 
entropy averaging schemes usually give similar results. 
The entropy average was not used in the present work. 


AREA AVERAGING SCHEME (AAS) 

The AAS uses conservation of mass (4), and area- 
weighted integrations of three quantities that might be 
measured directly by a slow- response probe behind a 
rotor, static pressure p, total pressure po, and total 
temperature to . 
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p = 7. 

/ pdy 
Jo 

(16) 
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Po = y 

f Pody 
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to = T 

/ tody 
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Primitive variables are found in a straightforward 
manner, except that the flow angle is found using the 
MAS. The AAS does not conserve mass, momentum, or 
energy. It results in the lowest efficiency and highest 
loss of the three schemes in the present case and is 
generally the least useful of the three. 


RESULTS 

Figure 1 shows relative Mach number contours from 
the unsteady solution at one instant in time. The inlet 
Mach number of about 1.2 produces a bow shock that 
stands ahead of the rotor blade and extends into the 
passage as an oblique shock. Further downstream, a 
strong normal shock stands between the rotor blades. 
Because of the large interblade- row spacing and the 
transonic flow within the rotor, the stator has little 
effect on the rotor, and the two rotor passage flows are 
nearly identical. The rotor wakes do have a large effect 
on the stator, however, and the three stator passage 
flows are quite different. Relative Mach numbers in the 
stator are transonic although absolute Mach numbers 


i 


* 


4 


Figure 1. — Unsteady relative Mach number contours. 



Figure 2. — Unsteady entropy contours. 


Figure 3. — Experimental contours of turbulent kinetic 
energy in stator (top) and computed entropy contours 
(bottom). 



Figure 4. — Velocity triangles for rotor wake and core 




Figure 5. — Unsteady static pressure contours. 


Figure 6. — Unsteady lift history for rotor (top) and 
stator (bottom). 


(not shown) are between 0.5 and 0.8. The rotor wakes 
are chopped by the stator blades and the remnants can 
be seen vaguely in the passage. 

Entropy contours in figure 2 show the rotor wakes 
more clearly (the shocks produce less entropy than the 
boundary layers and are not seen at these contour lev- 
els). The high-entropy fluid in the center of the rotor 
wakes migrates towards the pressure side of the stator 
blades. 

The enlargement of the stator region in figure 3 com- 
pares the computed entropy contours with contours 
of turbulent kinetic energy measured experimentally 
using laser anemometryby by Hathaway, et al. [16]. 
Although different quantities are plotted, both give a 
good indication of the wake positions and agree quali- 
tatively. 

The migration of the rotor wakes was originally ex- 
plained by Kerrebrock and Mikolajczak [19]. The ex- 
planation is summarized briefly below. Consider the 
rotor exit velocity triangles shown in figure 4. Relative 
velocities are denoted by W and absolute velocities by 
V. In the rotor frame of reference the rotor core ve- 
locity Wcore and the rotor wake velocity are 

roughly parallel to the rotor blade pressure surface. 
Rotor velocities can be converted to the stator frame 
of reference by adding the wheel speed, V = W + rfi. 


Ideally, Vcore points towards the stator at an appro- 
priate incidence angle, leaving tU wa t e with a lower ve- 
locity and a higher incidence angle than Vcore* The 
difference between U wa fc e and Ucore is called the slip 
velocity . This slip velocity carries the rotor wakes to- 
wards the pressure surface of the stator blades. The 
migration of the rotor wakes cannot be calculated with 
an averaging-plane approach but may possibly be cal- 
culated using Adamczyk’s average passage model [2] 
with additional modeling. 

Unsteady effects on the pressure field are shown in 
figures 5-9. Static pressure contours are shown in fig- 
ure 5. Again the rotor passages are nearly identical 
but there are large differences between the three stator 
passages. The pressure surfaces of the stator blades 
show large differences in loading. The suction surfaces 
have at least three small regions of low pressure spaced 
regularly along the blades. These low pressure regions 
are not clearly associated with the rotor wakes. 

Figure 6 shows unsteady lift coefficients for the two 
rotor blades and the three stator blades during one 
blade passing interval. The rotor lift varies by less 
than 1 percent, with three oscillations per blade pass- 
ing interval due to the three stator blades. The stator 
lift varies by nearly 30 percent, with two oscillations 
per blade passing interval due to the two rotor blades. 
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Figure 7. — Rotor pressure distributions at maximum 
and minimum lift, and time average of unsteady dis- 
tribution. 



0.0 0.2 0.4 0.6 0.8 1.0 

FRACTION OF STATOR CHORD 

Figure 8. — Stator pressure distributions at maximum 
and minimum lift, and time average of unsteady dis- 
tribution. 

Figure 7 compares the rotor blade pressure distribu- 
tions at maximum and minimum lift to a time-averaged 
distribution made by integrating the unsteady pressure 
at each point over time. The three curves are virtually 
identical. A similar plot is shown for the stator in fig- 
ure 8. Here, periodic fluctuations about the mean can 
be seen on both surfaces. On the pressure surface, the 
fluctuations have a long wavelength. On the suction 
surface, three distinct minima corresponding to the low 
pressure regions in the contour plots (figure 5) can be 
seen. 

* An interesting discussion of unsteady stator pres- 
sures calculated for a compressor stage can be found 
in the paper by Valkov and Tan [20]. They proposed 

* that the pressure surface fluctuations are caused by un- 
steady vortices formed on either side of the rotor wakes 
as they impinge on the stator pressure surface. They 
also proposed that the suction surface fluctuations are 
caused by small vortices formed as the rotor wakes 
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FRACTION OF STATOR CHORD 


Figure 9. — Comparison of steady and time-averaged 
unsteady stator pressure distributions (all averaging 
schemes). 



Figure 10. — Comparison of steady and time-averaged 
unsteady Mach number distributions across stator 
wake (energy averaging scheme). 


pass over the stator leading edges at a high relative 
incidence. These small vortices convect with the slow 
boundary layer flow on the suction surface and thus 
do not correlate with the wake locations. Although 
neither of these phenomena has been observed directly 
in the present work, the explanation is plausible. The 
large pressure fluctuations do seem to be a direct re- 
sult of viscous interaction between the rotor wakes and 
the stator boundary layers and could be a significant 
source of noise in turbomachinery. They are strictly an 
unsteady phenomena that cannot be predicted with a 
steady analysis. 

Figure 9 compares the time-averaged unsteady sta- 
tor pressure distribution with steady averaging-plane 
results made using each of the three averaging schemes. 
Despite the large unsteady pressure fluctuations in the 
stator, the time- average stator pressure distribution 
agrees almost exactly with the pressure distribution 
predicted by the steady analysis, regardless of the av- 
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ADIABATIC EFFICIENCY 


0.00 



adiabatic 

Averaging Scheme | 

efficiency 

energy 

momentum 

area 

unsteady 

0.88032 

0.87696 

0.87447 

steady 

0.87858 

0.86964 

0.85904 

difference 

0.00174 

0.00732 

0.01543 


Table 1 : Comparison of stage adiabatic efficiencies. 
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FRACTION OF BLADE PASSING CYCLE 


stator loss 
coefficient 

Averaging Scheme 

energy 

momentum 

area 

unsteady 

0.02616 

0.03236 

0.03582 

steady 

0.02809 

0.03042 

0.03256 

difference 

-0.00193 

0.00194 

0.00326 


Figure 11. — Stage adiabatic efficiency (energy averag- Table 2: Comparison of stator loss coefficients, 
ing scheme.) 


eraging scheme used. It is both surprising and encour- 
aging that the average pressure distribution can be pre- 
dicted so closely with a steady analysis. 

All of the considerations raised earlier regarding spa- 
tial averaging of a flow field also apply to temporal av- 
eraging. Again, nearly any quantity can be integrated 
in time with varying degrees of physical significance. 
In figures 7-9 for example, surface pressures were in- 
tegrated in time to produce average pressure distribu- 
tions. For the remaining results, unsteady stator exit 
conditions were first spatially averaged using each of 
the three schemes. Each scheme gives a complete but 
different vector of the average conservation variables q 
(eq. 1) at each point in time. To integrate these quan- 
tities in time, the integrands used for the EAS were 
constructed from q and integrated in time regardless 
of the averaging scheme used to generate q originally. 
It is by no means clear that this is the best way to time 
average unsteady solutions. 

Figure 10 compares the time averaged unsteady 
Mach number profile at the stator exit with the steady 
averaging-plane result made using the EAS. The pro- 
files agree closely across the wake but not across the 
core. The average unsteady results show high Mach 
numbers near the suction side of the blade and low 
Mach numbers near the pressure side resulting from the 
migration of the incoming rotor wakes. The steady pro- 
file is roughly through the middle of the time-averaged 
unsteady profile. The irregular shape of the unsteady 
stage exit profile cannot be modeled with an averaging- 
plane approach. 

The stage adiabatic efficiency was computed using: 


- (Poa/Poi)^ ~ 1 

(< 03 /< 0 l ) — 1 


(19) 


of blade passing interval (time) in figure 11. The un- 
steady efficiency varies by about ±1 point over a cycle. 

The time average of the unsteady solution is also 
shown in figure 11 and is compared with the efficiency 
predicted using the steady EAS calculation. The two 
results agree within 0.02 points in efficiency. 

Table 1 compares the stage adiabatic efficiencies cal- 
culated using each of the three averaging schemes. The 
EAS gives the highest efficiency and the best agreement 
between the steady and unsteady results whereas the 
A AS gives the lowest efficiency and the worst agree- 
ment. 

The stator loss coefficient was also calculated from 
the steady and unsteady calculations using the three 
averaging schemes, and the results are presented in ta- 
ble 2. The loss coefficient is defined by 

P02 ~ P03 

IjJ — 

P02 “ P2 

where 02 denotes the stator inlet. The stator loss coef- 
ficient is difficult to compute from the unsteady calcu- 
lation because the losses are referenced to interblade- 
row quantities pi and po 2 which are themselves un- 
steady. Here the EAS was used to define the stator in- 
let state 2 , then the three averaging schemes were used 
to compute the loss at 3 from both the unsteady and 
steady solutions. The EAS gives the lowest stator loss 
and the A AS gives the highest. The steady losses are 
slightly higher than the time-averaged unsteady losses 
using the EAS but are lower than the unsteady losses 
using the MAS or AAS. The EAS and MAS both give 
fairly close agreement between the unsteady and steady 
calculations. There is a somewhat larger difference be- 
tween the unsteady and steady losses calculated with 
the AAS. 


where ()i denotes the rotor inlet where reference con- 
ditions pox and toi are constant, and Q 3 denotes the 
stator exit. The stage efficiency is plotted as a function 


SUMMARY 

Three schemes for averaging computed flowfields 
in turbomachinery, the momentum (MAS), energy 
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(EAS), and area (A AS) averaging schemes, were de- 
scribed and compared. The MAS conserves global 
mass, momentum, and energy. For 2-D adiabatic flow 
it is equivalent to a mixed-out average, giving the flow 
properties after the wakes have mixed out far down- 
stream. It provides the correct available thrust but 
tends to provide low efficiencies and high losses. The 
EAS replaces conservation of axial momentum with an 
integration of the ideal total enthalpy. It is intended to 
give a good measure of the local total pressure and total 
temperature, quantities needed in the calculation of ef- 
ficiency and loss. It generally gives the best comparison 
with experimental data. The AAS integrates quanti- 
ties often measured experimentally: static presure, to- 
tal pressure, and total temperature. It has little other 
physical significance and gives the lowest efficiency and 
highest loss of the three schemes. 

The three averaging schemes were used as the basis 
of a steady averaging-plane analysis of a transonic fan. 
The rotor and stator were analyzed independently with 
average rotor exit conditions specified at the upstream 
boundary of the stator. The averaging-plane results 
were compared to an unsteady rotor-stator interaction 
solution which did not require averaging between the 
blade rows. 

Several interesting physical phenomena were noted 
in the unsteady solution. Large pressure fluctuations 
which could contribute to noise generation were seen 
on the stator blades. Despite the large pressure fluctu- 
ations the averaging-plane analysis matched the time- 
averaged pressure distribution almost exactly, regard- 
less of the averaging scheme used. The rotor wakes mi- 
grated towards the stator pressure surface. Although 
the effects of wake migration cannot be modeled with 
a averaging- plane analysis, the overall stage efficiency 
and stator loss coefficient were predicted fairly well 
with the averaging- plane analysis based on the EAS. 

Several unresolved issues remain concerning averag- 
ing schemes. First, many other schemes can be devised, 
such as mass- averaging total pressure for comparison 
to experimental data, or averaging entropy to estimate 
losses. These schemes must be evaluated for their par- 
ticular applications. Second, the fan stage analyzed 
has a large interblade spacing, so the results might not 
hold for a stage with stronger acoustic coupling. Third, 
the schemes can be extended to 3-D by addition of a 
spanwise integral, but it is not immediately clear what 
quantity to integrate. Finally, it is also unclear how 
to average unsteady solutions for the best comparison 
with steady solutions or with low frequency response 
probe data. 
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